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Abstract 

We conjecture that a non-flat D-real-dimensional compact Calabi-Yau manifold, 
such as a quintic hypersurface with D = 6, or a K3 manifold with D = 4, has locally 
length minimizing closed geodesies, and that the number of these with length less than 
L grows asymptotically as . We also outline the physical arguments behind this 
conjecture, which involve the claim that all states in a nonlinear sigma model can be 
identified as "momentum" and "winding" states in the large volume limit. 
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1 Introduction 



The two-dimensional nonlinear sigma model (NLSM) is a central topic in string theory, 
in statistical mechanics, and in math-physics interface topics such as mirror symmetry. It 
is a quantum field theory of maps from a two-dimensional (2d) Riemann surface into a 
Z)-dimensional Riemannian manifold M, a sort of quantized version of the harmonic map 
problem. Although all existing treatments of its perturbation and renormalization theory 
involve choosing coordinates on M, the physical results are covariant under diffeomorphisms 
PP, so the NLSM provides a direct contact between quantum field theory and geometry. One 
can supersymmetrize the NLSM and find even richer connections with geometry. 

In a sense which is still not well understood, the NLSM defines a generalization of Rieman- 
nian geometry which is often called "stringy geometry." While there are many interesting 
results in this subject, surveyed in [2], this paper will actually be more about the "large vol- 
ume" or "large structure" limit in which the NLSM can be understood using conventional 
geometry, and it will try to make statements about conventional geometry based on the 
properties of the NLSM. 

Our main conjectures are stated in the abstract and will be restated below. The first is 

Conjecture 1 A compact Calahi-Yau manifold has nontrivial closed geodesies which are 
local minima of the length functional. 

More generally, this should apply to any compact manifold for which the NLSM leads to a 
superconformal field theory (with H = 0, see section 2). The physics argument is simply that 
the NLSM must have stable states corresponding to strings winding about geodesies, but 
if a geodesic is not locally length minimizing, the corresponding winding string will not be 
stable. For this reason we will often refer to locally length minimizing geodesies as "stable" 
geodesies. 

From a physics point of view this claim may seem reasonable and unsurprising, but it 
has never been shown mathematically and there are geometric considerations that make it 
somewhat more surprising. While it has been shown that all compact Riemannian mani- 
folds have closed geodesies [5], in general such geodesies are not local minima of the length 
functional. For example, the sphere has no such geodesies; all are unstable at second order. 
More generally, as we review below, the second variation of the length is the negative of a 
component of the Riemann tensor, so positive curvature is an obstruction to stability. Now, 
since the Ricci tensor is a partial average over the Riemann tensor, and a Calabi-Yau mani- 
fold is Ricci fiat, at every point there is some two-plane with positive curvature along which 
(by analogy to the sphere) one might be able to vary the geodesic and lower its length. 

In the simplest example, the Eguchi-Hanson space, one can check that there are no 
stable closed geodesies, as we do in section 4. The only candidate is the geodesic winding 
the exceptional cycle (minimal volume nontrivial two-sphere), but it is not stable. Of course, 
the Eguchi-Hanson space is not compact, so it is not a counterexample to our claim. The 
simplest Calabi-Yau on which we expect to find stable closed geodesies is the resolution of 
the orbifold (M^ x S^)/Z2, which in a sense is two Eguchi-Hanson spaces glued together [HIS]. 
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This orbifold CFT has winding states which are the Z2-invariant projection of the winding 
states on x S^. By continuity under varying moduh, similar winding states must be 
present after a small deformation; these are the stable closed geodesies of our conjecture. 
While we do not show this explicitly, we will check that there are stable noncompact geodesies 
(escaping to infinity) on the Eguchi-Hanson space which might be glued together to produce 
these stable closed geodesies. Similarly, the K3 obtained by small deformation of T^/Z2 will 
have stable closed geodesies corresponding to the winding states on this orbifold, and so 
forth. 

The main physical point we need to justify our conjectures in general is to show that 
these NLSMs in fact have stable states corresponding to strings winding about geodesies. 
Our primary argument will be based on modular invariance of the partition function on the 
torus. This is an invariance under exchanging the A- and B-cycles of the torus, which in a 
sense relates the spectrum of the Laplacian on M, called the "momentum states," to the set 
of stable closed geodesies, the winding states. Showing this for the torus T'^ is a standard 
calculation, which we review in section 2. This argument can also be used for orbifolds of 
the torus, justifying the claim for small deformations of orbifolds. 

We will make some steps towards a general argument for the same claim, that modular 
invariance relates the spectrum of the Laplacian on M to the set of stable closed geodesies, in 
the large volume limit. Let us now recall some facts about this limit. In quantum mechanics, 
the quantization parameter (Planck's constant h) has (dimensional) units involving both 
length and time. It can be thought of as controlling an uncertainty relation between position 
and momentum, or energy and time; two related but distinct conjugate quantities. Analogous 
but different relations apply to most quantum field theories. However, as we will review in 
section 2, because of conformal invariance in two dimensions, the quantization parameter in 
the 2d NLSM has units of squared length on M, and controls an uncertainty relation between 
position and position. It is usually denoted as a' or in the string theory literature, and 
it determines the length scale Is = y/a' on M at which quantum fluctuations of the two- 
dimensional surface are important. 

The limit a' — )■ 0, or equivalently a limit in which we fix a' and scale up the metric on 
M by an overall constant, is the large volume limit. While it has many features in common 
with the semiclassical limit h ^ 0, there are some differences of interpretation which will 
become important below. A good starting point for explaining this is to describe the state 
space and Hamiltonian of the NLSM and compare it to the state space and Hamiltonian 
of the quantum mechanics (QM) of a particle moving on M. Recall that in most quantum 
theories, the state space is an infinite-dimensional Hilbert space "H, and the Hamiltonian H is 
an unbounded operator but with bounded below spectrum. Its spectrum and the "partition 
function" 

ZiP)=TTnexp-f5H (1) 

is a basic physical observable. In QM, the Hilbert space is the space of functions on 
M, corresponding to the quantization of a classical particle whose state is a choice of point 
on M and a conjugate momentum. The Hamiltonian is the sum of the Laplacian A on 
M multiplied by ff' and a multiplication operator by a function V on M (the potential). 
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The asymptotics of its spectral density and thus the partition function are governed by the 
Laplacian and are given by Weyl's theorem. Let us set the potential V = 0, then we have 



ZqmW) = Tr ^ exp -/3A ~^^o ^Vol (M) (1 + . (2) 

The Hilbert space of the NLSM corresponds to the quantization of a string (loop) on 
M, and should be some space of functions on a loop space on M. This type of definition 
has been worked out mathematically for M fiat (or parallelizable as in the WZW models), 
but due to the difficulties of renormalization this has never been done for other M. What 
can be done is the analysis in the large volume limit and the development of a renormalized 
perturbation theory in a'. We will review the NLSM spectrum and partition function in 
the large volume limit using physics techniques below. In the simplest case in which the 
geodesies are isolated, we find 

ZNLSKlij ) = {^Qm{(^ = Ta) + Zwinding{(3 = t/ tt')^ (3) 

Zujindingi/3) = ^exp-/3L(7)^ (4) 

7 

Here 77 (r) is the Dedekind eta-function or "classical partition function," and the sum in 
Zwinding is taken over all closed locally length minimizing geodesies 7 on M, while L{'-f) is 
the length of a geodesic. Thus the NLSM Hilbert space has two components. The first, often 
called the space of "momentum states," is the tensor product of two factors: the quantum 
mechanical Hilbert space of a particle moving on M, corresponding to a wave function of 
the center of mass of the string, and a Hilbert space of "oscillator states" which physically 
correspond to small fiuctuations of the string. 

The second component, the "winding states," is the tensor product of a Hilbert space 
with a basis vector for each closed stable geodesic 7, with a similar space of oscillator 
states. This component is not present in the quantum mechanics of a particle; its origin is 
what one would expect intuitively; a loop can wind about (embed into) a geodesic 7 to give 
a physical state. To correspond to a state, the geodesic must be a local minimum of the 
length functional; otherwise the state would be "unstable" and decay into a loop of lower 
energy (length). 

Of course, on the torus, geodesies come in families, and a priori there might be families 
of geodesies on a nontrivial Calabi-Yau manifold. In the physics, a family of geodesies will 
contribute a function to Zwinding, obtained by the "collective coordinate prescription." This 
amounts to finding the moduli space of the family, call this Ai, and quantizing the moduli 
space, changing Eq. @ to 



Zwinding{l3) = ^exp{ - Ttt'A^} (5) 

7 

where A_m is a Laplacian on the moduli space. 
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A similar computation and picture apply to the supersymmetric NLSM. The analog 
of Zqm is a supersymmetric QM partition function, which counts eigenfunctions of the p- 
form Laplacians for all p. The other terms Zwinding and l/|?7p'^ have analogs which will be 
discussed in section WA\ 

This picture is only known to apply in the large volume limit, as the existing justifica- 
tions involve perturbative expansions in a' which are believed to be asymptotic with zero 
radius of convergence. As one moves away from this limit, i.e., considers manifolds M with 
geometry on scales Ig or shorter, not much can be proven, but there are many physics claims 
for which a great deal of evidence has been assembled. The most basic of these is that 
supersymmetric NLSM's exist for certain non-flat manifolds M: the Calabi-Yau manifolds 
with SU{n) holonomy, and hyperkahler manifolds. It is also believed that supersymmetric 
NLSM's exist for manifolds M with the special holonomy groups G2 and Spin{7). 

The argument that NLSM's exist with Calabi-Yau target spaces, while not mathemati- 
cally rigorous, is extremely compelling, and has three parts [2]. The first part is that these 
NLSMs have unobstructed deformations, so that they come in moduli spaces of computable 
dimension. The second part is that there are supersymmetric conformal field theories, the 
Gepner models, which can be exactly constructued using algebraic techniques, and which 
can be compellingly argued to correspond to points in the moduli spaces of Calabi-Yau sigma 
models. Finally, the moduli spaces themselves can be explicitly determined using algebraic 
geometry and mirror symmetry. Thus, there exist families of NLSMs with both a large 
volume limit, and other "stringy" NLSMs which are not described by this limit. 

The properties of these "stringy" NLSM's are a primary question of the still-nascent 
theory of "stringy geometry," nascent because it has no general definitions or techniques 
at this point. A good deal of progress has been made on questions which can be answered 
in the topologically twisted NLSM, leading to the many results of mirror symmetry, and 
connections to mathematics such as quantum cohomology and derived algebraic geometry. 
But much less progress has been made on other questions, such as the general nature of the 
state space and partition function Eq. ([3]). 

Of course, finding the exact spectrum of the Laplacian or equivalently computing Eq. ([2]) 
exactly for general Riemannian manifolds is already an intractable problem, and it is not 
clear why computing Eq. should be any easier. Although from this point of view it is 
interesting that explicit expressions for Eq. ([3]) are known for Gepner models, our discussion 
here will not make use of this, but rather focus on qualitative properties. The main idea 
we will use is an analogy with the semiclassical trace formula [714T0]. In general terms, 
a trace formula relates the spectrum of the Laplacian on a Riemannian manifold M, to 
the lengths and other properties of closed geodesies on that manifold. Intuitively, such a 
relation will arise by taking the semiclassical limit of the functional integral over closed 
paths in M. This computes the trace of the heat kernel, which determines the spectrum, 
while the semiclassical limit is dominated by classical solutions, the geodesies. Typically, 
this relation is only asymptotic in h, except for special cases such as tori and homogeneous 
spaces. For the torus, the relation can be verified analytically using Poisson resummation, 
while for quotients of hyperbolic space it is the Selberg trace formula. 
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Looking at Eq. ([3]), it involves the same two ingredients, a sum over eigenvalues of the 
Laplacian, and a sum over closed geodesies. And, at least on a heuristic level, it is easy to 
relate it to the trace formula. Consider a semiclassical treatment of the genus one partition 
function; now the classical solutions are (by definition) the harmonic maps from the torus to 
M. This includes "world-sheet instantons," but it also includes simpler solutions in which 
(say) the r-cycle of the torus maps into a closed geodesic in M, and with no dependence on a. 
By the intuition which led to the trace formula, the sum over these solutions should be related 
to a sum over the spectrum of the Laplacian. In the usual discussion these are "momentum 
states," created by vertex operators which are local in M. Conversely, configurations with 
dependence on a but not r are winding states, and the genus one partition function also 
contains a sum over these. Now, modular invariance (or the "S'-transformation" r — — l/r) 
relates these two sums. Thus, the trace formula is part of the explanation of modular 
invariance in the 2d NLSM. 

By analyzing the modular invariance relation between the two terms in Eq. ([3]), we can 
understand the asymptotics of the number of geodesies as a function of their length. Let us 
first suppose for simplicity that the geodesies are isolated, then we will argue in section 14.51 
that the symmetry between the two terms of Eq. ([3]) will require Zqm and Zwinding to have 
the same asymptotics as r — )■ 0. This will imply that 

Conjecture 2' On a compact Calahi-Yau manifold of real dimension D, assuming geodesies 
are isolated, the number of nontrivial closed geodesies which are local minima of the length 
functional of length less than L, grows asymptotically as . 

To understand the case in which the geodesies are not isolated, one needs to know more 
about their moduli spaces. While we do not have much to say about the general case, the 
fact that Conjecture 2' holds in the torus and deformed orbifold examples, combined with 
the idea that the asymptotics of the individual terms in Eq. (jS]) and Eq. ([5]) cannot change 
under deformation to a nearby conformal field theory, strongly suggests that we do not need 
the additional hypothesis, so we make 

Conjecture 2 On a compact Calahi-Yau manifold of real dimension D, the number of non- 
trivial closed geodesies which are local minima of the length functional of length less than L, 
grows asymptotically as . 

We hope that further development of these ideas will allow making this argument more 
compelling, and perhaps determine whether geodesies on a general Calabi-Yau are isolated 
or not. 

To conclude this introduction, let us mention some loosely related work. In [11] it was 
suggested that NLSM flows between target spaces which are higher genus Riemann surfaces, 
and Liouville-type theories, could be defined by using the Selberg trace formula to compute 
sums over winding states. In [12], a general relation was proposed between modifications 
of the contour of functional integration in a d-dimensional QFT, and boundary problems 
in a topologically twisted d + 1-dimensional QFT. It is tempting to imagine that relations 
between d = 1-dimensional trace formulas and d = 2 superconformal field theory could be 
understood in these terms. 
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2 Winding states in sigma models 



We begin by discussing the bosonic nonlinear sigma model with fields : S — t- which 
are local coordinates on M, and the action 



^ = -i_ J d?xy/hh''^da,X^'dpX^G^^{X) + 1 j X*{B). 



(6) 



Here h^ii is the worldsheet metric and G^i^ the metric on M. The indices are = a, t; 
and /i, V are space-time coordinate labels. 

The i?-field is a local two- form on M and X* is its puUback to S. In general B need not 
be globally defined - its contribution to the action is only defined up to shifts of 27r, and in 
the supersymmetric case anomalies can enter. Nevertheless the three-form H = dB will be 
globally defined. It will be zero in the NLSMs we discuss, but in this section we work out a 
few results for 7^ 0. 

In a semiclassical treatment of the functional integral, we sum over integrals defined by 
expanding around solutions of the classical equation of motion. This is the harmonic map 
equation, generalized by i/-fiux, 

= d^iVhN'^dpX^) + \fhK'PTl^^{X)d^X''dpX'' + e"^G^'^(X)/7,A<x(X)9,X^9^X'^. (7) 

where r^_^ is the Christoffel symbol for the metric G^j^p. It is independent of the two- 
dimensional conformal factor. 

We now take h"^ = 6°^^, and consider the special case in which the fields only depend 
on one worldsheet coordinate, say a. In this case the if-fiux drops out, and these classical 
solutions are closed geodesies on M. 

Let 0" G [0, L) = /, and take 

-f-.I-^M (8) 

to be a closed geodesic on M, satisfying Eq. ([7]) and 7(0" + L) = 7(0"). The equation Eq. ([7]) 
implies that 

= ^|7'(a)r (9) 

where 

|y(a)P = G,.(7(a))^^, (10) 

so up to a factor, the geodesic is parameterized by arclength. We now choose 

l = |7'(a)p (11) 

so that the geodesic is parameterized by arclength, and L is its length. 

Since conformal transformations necessarily mix a and r, a nontrivial geodesic breaks 
conformal symmetry. 
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2.1 Brief review of the torus target space 

Let us now review the well known case of M = with constant metric G^^ and B^^. We 
take the coordinates X'* to range over the unit hypercube [0, 1)-^. Thus the closed geodesies 
are 

X'' = x^ + w^"- (12) 
Li 

with xl e [0, 1)^ and wf" e Z^. 

To get the states of the quantum theory, the zero modes Xq must be quantized, leading 
to discrete momenta valued in the dual lattice 1,^. 

The energy of a winding state is simply Eq = J da^^. Taking into account the kinetic 
energy the torus partition function is 

Z{t) = ^ exp 1 7rr2- {-JY~ + ^^1^ ' A^) + 2i7rriW ■ /i| (13) 

In this form, the partition function has a natural interpretation as a trace over Hilbert 
space. This can be seen by recalling q = e^'^*'^ and 

Z(r)= g^^^+^^^V^^^-^^)' (14) 

Poisson resummation formula makes this explicitly modular invariant 

(27r)^ ^ _ f ttLV-™" 



= L I \ (15) 



Adding the oscillators further contribute a power of the Dedekind function |?7(r)|~^^. As 
in well known, the above sum is invariant under both shift (T : r — )■ r + 1) and inversion 
(5* : r — )■ — 1/r) which generate SL{2, Z). 

The torus is special in that each winding state also carries momentum. This type of 
degeneracy does not arise for an isolated closed geodesic on general curved spaces, and in 
such a case the winding or momentum quantum numbers do not occur simultaneously for a 
sector of the Hilbert space. Modular invariance does not come as a trivial consequence of 
summing over orbits of the SL{2, Z) group action. 

2.2 Expansion to second order around a harmonic map 

To generalize this to curved M, we want to expand the action Eq. around a solution Xq, 
schematically 

X(a,r)=Xo(a,r)+e(^,r). (16) 

We are free to define the coordinates C,{<y, t) in any way we wish so as to simplify the 
expansion. 
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In the usual covariant treatment of the sigma model, one expands around X fixed to a 
point p, and takes ^ to be Riemann normal coordinates (RNC) around p. These are defined 
by considering the geodesic fiow (or exponential map) starting at p; the coordinate of a point 
q is the initial velocity ^ of a geodesic which reaches q at time t = 1. Thus, geometrically, C, 
is a tangent vector at p, in other words C, G TMp. 

A natural generalization is to take C,{cr, t) to be a RNC around the point Xo(cr, r). Thus, 
C.^icr, r) is a tangent vector to the space of maps X : S — )■ M, which can be regarded as a 
lift of the map ^ : S — )■ TM satisfying vr^ = Xq. 

To quantize to one loop order, it suffices to expand the action to second order in the 
fiuctuations ^. There is a simple covariant result [I3] for the expansion to this order around 
an arbitrary solution X of Eq. ([7]) with H = 0. It is 

S = Jd'a Vhh'^^ [G^,{X{a, r))^^^^^^ - Rx^.pAX{a, T))d^X%XPeC] ■ (17) 

In the simple case of constant X{a,T), the curvature term drops out, leaving the leading 
non-interacting term in the usual a' expansion. The curvature term is new and arises at the 
same order {i.e., it is independent of a') in the process of covariant izing the second variation. 
Note that both terms are independent of the two-dimensional conformal factor. 

A simple geometric way to compute this is to consider a family of maps depending on 
two extra parameters {u,v), 

F : X S ^ M, (18) 

of the form 

F = X + <i + <2, (19) 

and then take the first variation with respect to each of the new parameters. Varying 
the metric will give connection and curvature terms, in a very analogous manner to the 
fermion connection and curvature couplings generated in the superfield formalism by the 
^-dependence of the metric. Taking conformal gauge, we have 

6J,S = j d^a 2G^,(X)9e^9r + ^G,,,,x{X)ededX' + G^,,x,{X)i''i''dX^dXr (20) 

We then need to rewrite the second term in terms of (^)^ and (9^)^. Taking the symmetric 
and antisymmetric combinations, we have 

4^^^^^^ = 2d{i^i^) + 2 [i^di^" - i^'di^) . (21) 

The symmetric term can be integrated by parts, to produce another G^^^xa term and a term 
with d'^X^ . This can be written using Eq. ([7]) T{dX f termQ The final result is 

- 2G,,,A/.(^)eV5^''5^' + 0{T) + 0{TH). (22) 

satisfies the classical equation of motion by assumption. 
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As for the antisymmetric term, if the final result is covariant, there is no antisymmetric 
tensor we can make out of R{dXY . Furthermore, if we grant that the final result is covariant, 
we can easily get it by using Riemann normal coordinates at X{a, r), in which G^^p^x = Tj^^ = 
0. We then have 

2 

G ^u,\a — '^Rfj.Xua (23) 

which combines with the 1 — (—2) above to give —R as expected. 

A similar computation can be done for the X*{B) term. The first term combines with the 
antisymmetric part of Eq. fl2T|) (after integrating by parts) to produce the expected Bd^d^ 
term. The last term combines with the symmetric part to produce d\Hfj_yfjdX^dX'^^^^" . 
This term is covariantized by the TH term coming from the equation of motion. In addition 
one finds a H^d^dX term. 

Thus the second order variation with H is the sum of Eq. (ITTl) and 

S = j (fa e"'^ [5^.(X(a, T))d^edpC + H^^^^d^X^^dpC + V ,H,xp{X{a, T))d^X^dpXPeC] ■ 

(24) 

We see from this that if Vi/ 7^ 0, it can contribute to the mass term, and thus the definition 
of stability of a geodesic will change. This is relevant for Wess-Zumino-Witten models, for 
example, where it allows for stable geodesies on group manifolds with positive curvature. 
We will assume H = from now on. 

This expansion becomes complicated at higher order. For X a geodesic solution, this can 
be simplified by using Fermi normal coordinates, as we discuss in Appendix ( lA.ip . 



2.3 NSR superstring 

The generalization to the NSR superstring is very similar as we are expanding around ip = 0, 
so we just use the standard fermion action, 

Sf = j fa {6-;d, + d^X'rQ i;^ + ^R,xu.Vi;^rr ■ (25) 

The curvature term is not relevant at one loop, and since F = along the geodesic in Fermi 
normal coordinates, the action becomes free. There is a nontrivial boundary condition 
determined by the holonomy of the geodesic. 

In the Appendix flA.2l) we show that the nontrivial geodesic also breaks worldsheet su- 
persymmetry. There is however no fermionic zero mode in the Neveu-Schwarz sectors and 
so these susy breaking geodesies will contribute as saddle points of the path integral. 

A remnant of the broken susy is that the longitudinal bosonic degree of freedom and 
the corresponding (real) fermionic fluctuation, preserves one supercharge that acts trivially 
on the transverse modes (see Appendix IA.2p . This remainingsuper charge clearly does not 



accommodate for example spectral flow of the original theorjo, so in such winding sectors, 
the Ramond and Neveu-Schwarz states are no longer connected. 



^There is no U(l) symmetry, at most 0(1) which is a sign change. 
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In other words, the partition function including contributions of these saddle points need 
not transform nicely under the spectral flow. This is no accident, the continuous change of 
moding by the holonomy of the winding string is a physical input and not a choice. Phrased 
in another way, the partition function is not topological. 

It is worth emphasizing that holonomy caused by Poincare map does not effect the 
supercurrents moding numbers, as will become more evident in section (13. 3p . As a result, as 
expected these (non-BPS) states does not contribute to for example the Witten index. 

2.4 Mass terms, positive and negative 

For a general metric and expanding around a general solution, even computing the propaga- 
tor for the quadratic action Eq. ( IT7|) is a difficult problem. It can be simplified somewhat in 
the case of a geodesic by the use of Fermi normal coordinates. These are defined by expand- 
ing the tangent vectors C, in terms of an orthonormal frame, defined by parallel transport 
along the geodesic. 

Their main properties which we will use in the following is that 

Gfiul-y — ^fiu'i r^j^l^ = 0. (26) 

One may ask is there a global obstruction to doing this for a closed geodesic? What about 
the overall rotation of the basis? For example, one may expect the overall holonomy to be 
trivial, as the closed geodesies we are interested in are topologically trivial. This however 
need not be the the holonomy characterizes the behavior of the cotangent bundle, 

as will be clear momentarily (see section (13. 3p ). 

Thus, in these coordinates, Eq. ( IT7|) has a canonical kinetic term, and describes massive 
scalar fields with a position-dependent mass. Since the action is independent of t, we can 
go to normal modes ^ oc exp(mt). The resulting equation of motion is a time-independent 
matrix Schrodinger equation, 

-^e^ + M,^(a)C = n2e^ (27) 

with 

Mna)^-R^,^M^))iy)\ir- (28) 

Here indices are raised and lowered with G^j, = 5^,^, so M^,y is a symmetric matrix. Note 
that 

M^^iyr = Vi. (29) 

so the longitudinal fluctuation is massless and decoupled from the other fields. 

With n = 0, this is essentially the geodesic deviation equation. If we start at ,^ = and 
slightly vary the initial velocity the qualitative behavior depends on the sign of M (and 
thus R). For M < 0, the solutions will be oscillatory, and nearby geodesies will stay nearby. 
This is the case of positive curvature, such as a sphere. On the other hand, for M > 0, ^ 
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will grow exponentially and nearby geodesies will diverge from each other. This is the case 
of negative curvature. 

When considered as a winding mode in string theory, the dynamics is a bit different. Here 
M is playing the role of a world-sheet mass term. For M > 0, i.e. negative curvature, the 
bosonic fields are massive, and conformal invariance is broken. For M < 0, i.e. positive 
curvature, the system appears to be unstable. This is the reverse of the previous discussion, 
which at first may seem a bit paradoxical. 

In fact there is no contradiction. In general terms, positive curvature allows a small 
variation of a closed geodesic to decrease its length, as is intuitively apparent for geodesies 
on the sphere. Conversely, negative curvature stabilizes the geodesies. Of course, on a Calabi- 
Yau manifold, the Riemann tensor will always have components of both signs (since the Ricci 
tensor is zero), and is furthermore not constant, so the discussion is more complicated. 

This observation is directly relevant for the Hamiltonian quantization of a winding state. 
If the matrix Schrodinger operator in Eq. fl27p has a negative eigenvalue, then the corre- 
sponding state should be unstable and is not in the spectrum. Later we will check this in 
solvable examples such as Eguchi-Hanson. It is possible that the time-dependent eigenvalue 
will cross zero multiple times, this means the instability is turned on for a fraction of the 
tim^. Then it would seem the corresponding winding state could still exist, however, see 
below. 

What is the order of the two terms? Because we are working in an orthonormal frame, 
the circumference of the a direction is L, and the normal modes in a will have rough energies 
VP? jL?. Components of the curvature in an orthonormal frame are roughly I/t^^^^, where 
Tf-urv is the curvature length. These are comparable in the simplest cases, but not in general. 

It is easy to find examples with L » Vcurv The simplest is to consider the n'th iterate 
of the geodesic, in other words the map defined by composing a — >■ na with 7. Although 
we are in a loose sense expanding around the same solution, the stability analysis changes 
because we allow ^ satisfying different (weaker) boundary conditions. The relative scaling 
suggests that if there is a region of the geodesic with negative M, then past some n the n'th 
iterate of the geodesic will be unstable. This is because the M energy will be proportional 
to n, while the level spacing will decrease as l/ra^. 

There is an additional aspect of the stability issue, which will be clearer after discussion 
of the next section. This involves the rotation of the normal coordinate around 7. Concrete 
example in section (14.41) shows that such contributions can overcome the negative mass term 
(from positive sectional curvature). All in all, the dynamics of the sigma model distinguishes 
our stability criterion from that of a geodesic in target space, and stability of a winding string 
is more intricate than the mass terms signs immediately indicate. 

•^For the specific example, Eguchi-Hanson space which we check in section (|4.4p . this possibihty is not 
reahzed. 
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3 Brief review of the trace formula 

The physical approach to the trace formulas is as a statement about the semiclassical limit 
of quantum mechanics [8]. Intuitively, they are based on the computation of a quantum 
mechanical partition function using a functional integral over closed classical orbits of a 
particle. This leads to a relation of the general form 

closed paths ' 

where S is the action of a closed orbit of time t, 1/ det is the one-loop approximation to the 
functional integral, and v is the Morse index of the orbit (the number of unstable directions). 
See for example chapter 17 of |H], or the reviews [TOlIT^ for precise formulas. 

A great deal of work by both physicists and mathematicians can be summarized in the 
general statement that formulas like Eq. ( 130|) are fairly well understood when treated as 
asymptotic expansions in but exact results (at finite K) are much harder to come by. The 
prototype is of course the Poisson resummation we used above in (fT5|) . while the most famous 
example (which initiated this field in mathematics) is the Selberg trace formula [7]. This 
applies to homogeneous spaces, such as a higher genus Riemann surface with a metric of 
constant negative curvature. Its standard proofs have little to do with the physics intuition 
and rely more on representation theory of groups. 

At the level of an asymptotic series in /I, then the physics argument according to which 
the functional integral can be treated by stationary phase, has (more or less) been made 
rigorous, as explained in [lOlIll], for the case of a point particle. We now give some more 
detail on the relevance (and distinction) of the trace formula in the context of NLSM, already 
alluded to in the introduction. 



3.1 Basic relation to modular invariance 

Suppose we have an asymptotic expression for the wave trace, 

/(t) = 5^e-**v/^, (31) 

n 

then we can get the heat trace from the integral 

Now, if we know that /(t) is analytic and falls off in the upper half plane, with poles at 
t = Ti with residue a^, we can evaluate this by contour integral to get 

E^"''" = 7^E"^^"'^^'''- (33) 
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Thus, the sum over momentum states Zqm appearing in Eq. ([3]), is related to the sum 
over winding states as 




y/T' 



r = 



1/r. 



(35) 



(34) 



Thus, taking r — )■ 1/r roughly exchanges the two terms. This relation r' = 1/r has then 
has the role of a modular transformation in the conformal field theory. The corresponding 
discussion regarding the nonlinear sigma model can be found in in section (14. 5p . 

3.2 The wave trace and Laplacian spectrum 

For our purpose, the following formula establishes a whole tower of (quasi-)eigenvalues of 
the Laplacian, for each periodic geodesic. This follows from studying the wave equations, 
the formula first appeared in the work of [T6] . 

Take the wave trace as defined in (!3T!) . expanding around a point of its singular support, 
the residue is expressed as a sum over iterations of a 'primitive' geodesic 



where A„ is an eigenvalue of the Laplacian operator on M, z/^ the Morse index of the geodesic 
7, T its period (i.e. length L^) and is the Poincare map represented on the cotangent 
bundle in a canonical normal basis. 

There are three different possibilities for the linearized symplectic transformation on 
the cotangent space, invariantly characterized by the Birkhoff quadratic form [15], the so 
called elliptic, real and complex hyperbolic. The eigenvalues of P^ follow the usual pattern 
of Sp{2n) matrices. For elliptic case one has complex conjugate pairs of phases, (e*^%e~*^') 
corresponding to rotations, we do not need the details of the other cases here. The crucial 
difference between the elliptic and hyperbolic cases is that the former leads to discrete 
Laplacian spectrum (or wave group spectrum as used in [15]), but not the latter. As we are 
interested in a stable geodesic on a compact space (hence discrete Laplacian spectrum), we 
focus on the elliptic case from now on. 

For such geodesies the P^ eigenvalues all have unit norm which can be represented by 
rotations on the position- momentum planes. Label the corresponding angles by di, then 
from the expression of g-«^/^* j^iay be extracted (A is second order) as in [T6] . 



where {D — 1) is the number of transverse (real) dimensions to 7 and n, Ui are indepen- 
dent integers. The role of this formula in relating momentum sector and winding sector of 




(36) 



27r |det(/-P^)|V2 




(37) 
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string states will be discussed in section (14. 5p . We mention that the approximate (quasi- 
)eigenmodes associated to these values can be constructed as a local solution, a 'Gaussian 
beam' which travels along the geodesic (for details see [T7]). 

In the case of momentum modes, there is a well known universal behavior for the asymp- 
totics. Weyl's formula of the Laplacian spectrum. Ignoring numerical factors, it states that 
the integrated density of states 

N{X<S)~^pU0(/-i') (38) 

Treating this as a phase space integral then determines the asymptotic of the Laplacian 
eigenvalues 

'^-{hk) '''' 



where k is integer. In terms of diameter then this is the expected {j^Y behavior inside a 
'box' of size L. As already explained this size may differ from the length of a general closed 
geodesic. We can compare this with (137|) . while the scaling agrees assuming the geodesic 
length being close to the manifold M diameter, there are certainly more parameters in 
(13 7p than in the universal Weyl formula. 

This contrast is even more evident when one realizes the Weyl formula is really universal, 
including the dropped numerical coefficient (which is the same as on M^). The leading term 
in the Weyl formula f l38|) does not depend on the point chosen on M, the metric on M 
or even M itself. The wave trace singularity on the other hand, depends on a lot of local 
information in the neighborhood of L^. 

There is in fact neither a contradiction nor a puzzle here, but a short mention of the 
relation between the two behavior which are both 'asymptotic' is perhaps clarifying, it is 
simply that the Weyl formula thought of in the right way is the wave-trace formula for 'zero- 
length' geodesies. We refer readers interested in the technical details to (chapter 8 of) [18]. 
For our subsequent physics discussions, we are mostly concerned with closed geodesies of 
non-zero length, and the Laplacian spectrum will be tacitly taken to be those arising from 
the union of the eigenvalues given by ( 137|) for the non-zero length closed geodesies. This is 
one of the 'improvements' of the stringy geometry versus the particly geometry. 



3.3 Poincare map and oscillator frequency 

The Poincare symplectic map as discussed in the wave trace formula is responsible also 
for inducing nontrivial boundary conditions for the (transverse) worldsheet fields z = 
1, 2, . . . , (D — 1), leading to shifted frequencies of the oscillator modes. 

To see this, we recall the geodesic may be considered an integrable Hamiltonian fiow, 
when we take the particle Hamiltonian to be 

n = \g,.{x)p^p^ (40) 
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which turns the 2nd order geodesic equations into the corresponding Hamiltonian equations. 
In Fermi normal coordinate along 7 this is simply the Laplacian A = l^fjrgfir- 

Due to this realization, the linearized Poincare map action on T*M preserves the sym- 
plectic form, being a canonical transformation. Making the assumption of an elliptic closed 
geodesic with no resonances @, we can find a (complex) basis in T*M which diagonalizes 
as phase rotations 

P^F, = e^'^Y,, = e-''% (41) 

where i = 1, 2, . . . , (D — 1). These are usually chosen to be Jacobi vectors along 7 which are 
evolved along the flow (see [19] )• In formal terms, these vectors span an isotropic subspace 
of the cotangent bundle (an A- type brane), and it is interesting to consider the open string 
version of our discussion with the annulus amplitude replacing the torus partition function. 

This basis will be different from the local parallel basis, and the conversion involves a 
unitary conjugation, it can be shown (using Wronskian of the Yi basis) [T7] that this leads to 
the so-called "Birkhoff normal form". Using our Fermi normal coordinate, this arises from 
the sectional curvatures (involving planes (Oi) and (Oj) say), which are in FNC the second 
derivatives of the longitudinal metric component in transverse coordinates q^q^j ■ 

For our winding string, there is a similar effect from the Poincare map on the periodicity 
condition which must be imposed on the winding states. Consider the local eigenbasis of 
the mass matrix M^j, which are the sectional curvatures for the tangent Oi-planes, i = 
1,2, . . . , [D — 1). To leading order, the holonomy around 7 contributes a uniform kinetic 
term due to the worldsheet action. Let the rotation angles be 6i, for the length geodesic, 
we find 

Equivalently we can see this from the canonical quantization of the worldsheet fields, taking 
the conjugate momenta to ^*(cr, t) as 



(a,t) = G(X),,9.e + ... (43) 



where . . . stand for fermionic terms vanishing identically in Fermi normal coordinate. 

The zero modes of (vTj, ^*) are simply the particle position-momenta pair and under the 
geodesic flow rotates by 



which are but the annihilation and creation operators and a, 



^Resonance happens when the frequencies of the harmonic oscillators in different transverse directions 
have integer linear relations, which leads to divergence due to small denominator. 

^Notice for discussion in the winding sector, we defined the conjugate momenta using a as time, instead 
oft. 
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We can introduce a formal infinite dimensional symplectic form for the fields , or we 
simply make all the oscillators get the same phase under pj^. As a result, the oscillator 
frequencies are shifted for the bosons to and the mode numbers shifted to The 

zero point energy for these winding sectors are then shifted as well, to 



]_ 1^ M 
24 ^ 47 W 



(45) 



which agrees with the estimate fH^ . The effects on fermionic fields are identical, and we 
emphasize that this is not just the curvature induced mass term, which is not present for 
fermions at the quadratic order of the action expanded in FNC. Also, as alluded to at the 
end of section (12. 3p . as the bosonic and fermionic fields shift in the same way in terms of 
oscillator modes, whether the supercurrents have zero modes or not are only determined by 
the sector the fermions belong to 'originally', i.e. Neveu-Schwarz or Ramond0. These shifts 
account for the various 6'-function contributions when we calculate the worldsheet one-loop 
determinant contribution to the torus partition function of the NLSM in the next section. 



4 One loop NLSM partition function 

In this section, we study the fluctuations around a stable closed geodesic on the target space 
M. Such fluctuations are the curved space counterpart of oscillator modes which build up 
the spectrum in each winding sector. We only consider the expansion of the NLSM action 
as in (12. 2 p to quadratic order, sufficient for the one loop approximation on the worldsheet. 
The contribution to the action by the classical geodesic solution 

^/<*^4vm^- (46) 

where (ni,n2) are the two winding numbers around 7(cr), i.e. the homotopy class of the 
harmonic map. In expanding around the winding sector, the sum then has a natural inter- 
pretation as a iteration of the coprime homotopy classes with gcd(ni,n2) = 1. Especially, 
for winding states the sum goes over (n, 0) only. 

The one- loop determinant comes from the action Eq. f fT7|) . which for h""^ = d'^^ and for 
X((T, r) = 7(0") becomes 

S = jd'a [G,M^m'dC - Rx,pMa)){i)\in'C] ■ (47) 

^Again, under Hamiltonian flow, the Hamiltonian equations dictate that + iir') — >■ + itt) where 
"s are the 'geodesic time' derivative and after integration around 7(cr) this is how the phase rotation is 
picked up. For the string case, we replace C by ^ which accounts also for the normalization change for 
string-oscillators . 

^The shift of moding is reminiscent of a spectral flow, which does not break supersymmetry. However, 
here susy is broken due to the explicit mass matrix for bosons and but not fermions. 
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The fermionic action is free. We have ehminated the hnear terms dXd^ as usual, employing 
equations of motion for X, however, we have nontrivial periodicity for the bosons and these 
linear terms contribute non-vanishing terms. This will be best appreciated in section (14. 5p 
where we recover them through modular invariance. 

The rest of the section is as follows. In section (14. ip we give an overview of the differ- 
ent contributions from the various bosonic and fermionic quantum fluctuations, taking into 
account nontrivial holonomy induced boundary conditions. In section (14. 2 p we explain the 
geometric significance of the action f H7|) for a closed geodesic. There discuss the zero modes 
of the quadratic action and we give an argument for the asymptotic energy of the oscillator 
modes in general curved space. 

We then make a further perturbative expansion of the quadratic action in the curvature 
tensor in section (14. 3p . allowing us to compute the leading contribution for a general Ricci- 
flat target space M. We also then comment on the subleading corrections in this expansion, 
relating them to sub-principal wave invariants introduced by [SU] [ID] in the context of 
spectral geometry. 

In section (14. 4 p we study an explicit (non-compact) example, the Eguchi-Hanson metric 
and discuss the stability of all geodesies on this space. In (14. 5 p we close the loop of ideas by 
showing that the trace formula fl37j) gives us a Zqm that agrees with Zuimding by performing 
explicitly the r — ?■ — l/r transform. 



4.1 Free and massive contributions 



We now explain the one loop determinant due to quantum fluctuations of various longitudinal 
and transverse modes. First the longitudinal bosonic mode contributes as a (non-compact) 
free boson at one loop order (worldsheet action) 



7 _B, longitudinal , 
^7,1— loop ' 



1 



r2\v{rW 



(48) 



This factor is modular invariant by itself and so is not very useful for inferring information 
about winding sector states. 

The fermions are also massless at one loop order, and we find the free fermion partition 
function with twisted boundary conditions (which we keep general) 



^7,l-loop(''") 



r/(r) 



1=1 



ai/2 
A/2 



(0;r) 



(49) 



where we write explicitly only the chiral sector. The explicit value of (aj,/3j) depends on 
periodicity induced by action of the Poincare operator along 7(0"). And ai represents a 
spatial twist (i.e. winding mode) while A is a gauging or temporal twist. The contribution 
of each (combined complex) fermion 



-e 



ai/2 

m 



(0;r) 



(50) 
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is known to transform as 



,i(r + l) = Z''„„(-l/r) 



(51) 



The phases under r — r + 1 constrains modular invariance in the sum over spin structures 
in a string theory partition function. In our case, the fermionic phase and the bosonic one 
cancel, as both have identical frequency shifts due to the Poincare map. 

It is perhaps amusing to consider behavior of the partition function when D = 4, letting 
the curvatures R^upx blow up along the geodesic neighborhood so that the transverse massive 
bosons decouple, then combining the remaining degrees of freedom we find the following 
result 1^ 



'7F,Bo / N 



6162 



(52) 



which bears a formal resemblance to the massive character^ often seen in computing K3 
elliptic genus [22] [23] . 

On the other hand, the full one-loop torus partition function for winding states including 
the transverse degrees of freedom (written for the case D = A) 



loop 



T 



E 

7 stable 
nGZ 



n 



=1,2 



(0;r) 



Hdet [{dl-e)l + M{a)] 



(53) 

The ^-function characteristics are determined both by the frame holonomy around 7 and 
the fermions fiat-space periodicity labelled by spin structure {si,Si), while most often we 
will restrict to NS sectors in our discussion. The infinite product over mode number A; is a 
reduction along the direction transverse to the geodesic direction on the worldsheet ('time' 
or t), the resulting 'index form' is often used in Riemannian geometry and will be the topic 
of next section. 



4.2 Zero modes and loop partitions 

This section is a small detour and explains some useful facts about Jacobi vector fields and 
the index form. The Jacobi vector fields are zero modes of the temporal reduced winding 
sector sigma model action, the index form. They give the moduli space for non-isolated 
geodesies mentioned in the generalized version of Claim (2). 

For an isolated stable geodesic, all fluctuations in the transverse directions are massive. 
In such cases, there are no smooth Jacobi vectors along 7. When 7(cr) belongs to a pa- 
rameterized family, it is useful to introduce the Jacobi vector field. It satisfies the following 

^Note the labels below in the 9i and 62 are only reflecting which holonomy condition they should be 
associated to, and are not the usual nomenclature of 0- functions as used in [21j . 

^More precisely massive characters of A/" = 4 supersymmetric su(2) current algebra at level 1 in the 
Neveu-Schwarz sector. 
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equation 

Vy Vy J(a) + R{y, J)y = (54) 

where in components the RHS is R''kij{l')''J''{l')i- We see that the Jacobi fields are the zero 
modes of the action. In fact its contribution to the action is a total derivative. 
Consider the following functional of 7(cr), usually called the index form 

/id = fda [G,,{^{am"e' - ^A.p,(7(^))(7')'(7')Te] (55) 

This is the reduction of our 2d quadratic sigma model action along 7. The Jacobi equation 
( l54l) immediately shows that Jacobi vectors minimizes the index form. The null space of the 
index form considered as a bilinear form on T^M is spanned precisely by the Jacobi vector 
fields. For a Jacobi field, applying (!54|) gives 

jda^ (G.M^))e^^ = {J, J')\o' = (56) 

The index form is positive definite for length minimizing geodesies. While if the geodesic 
7(0") contains conjugate points, the dimension of subspace in T^M on which the index form 
is positive definite jumps by an integer whenever it goes past a conjugate point (of a chosen 
origin) [25]. 

This integer is the number of distinct minimal geodesic between the conjugate pair. This 
defines an index, which is a monotonic function of the arc length parameter a. For a closed 
geodesic we then recover the conclusion from physical considerations in section (12. 4p . for the 
stability of the iteration of an unstable geodesic gets worse and worse as it wraps around 
more times. 

The Jacobi vectors are null, meaning it is orthogonal to all other vectors of TM with 
respect to the bilinear form /id. In other words, taking a linear combination of a non-zero 
mode (p{(j) and any Jacobi vector field J(o"), 

Cia) = ^\a) + J\a) (57) 

the action 

/id[e] = jda W'f - /?^,,,(7')^(7')>V^] = hM (58) 

since the cross terms vanish by Jacobi equation. This shows that such zero modes completely 
decouple from non-zero oscillator modes, and their quantization leads to the additional factor 
in ([5]) which measures motion on the moduli space of nearby geodesies. This aspect of our 
discussion is formal without an actual metric which supports non-isolated geodesies of such 
properties. 

From the perspective of Morse theory [25] , the one dimensional sigma model is the Hessian 
of the energy functional whose critical points among all closed loops are the closed geodesies. 
The loop space VLM is the configuration space of the quantum sigma model. This allows 
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us to deal with semiclassical modes geometrically, by thinking of them as a (transverse) 
standing wave dissecting the closed geodesic 7 into n strands. In the case of an isolated close 
geodesic, the above are the favored modes of motion (up to longitudinal reparametrization). 

As a result we find a n-piece partition of the geodesic 7(0") by a nearby loop Un, not 
necessarily geodesic itself. For a fixed geodesic 7 G fiM, beyond a certain mode number N 
the partition is fine enough that generically all the (string) strands are each length minimizing 
geodesic arcs. On M the length scale determining the minimal for which all n > 
partitions are fine enough is the injectivity radius, r^-^y It does depend on the partition 
points. 

Taking the minimal value of injectivity radius rinj along 7(0"), label it rinj^^ We can 
estimate the minimal partition number to be 

N,^^ (59) 

' mj 

For n > then each strand is within injectivity radii of all the partition n-tuple of points 
(pi,P2, ■ ■ ■ ,Pn) = 7(0") n Un involved. Hence the index form /id is positive definite on each 
interval (uj, Cj+i), moreover its value is bounded below by that of the geodesic arc ■yai,ai+i- 

Hence the growth of energy for the partition for a perturbation of 7 G QM to a loop Un 
nearby, is bounded below by a sum of n terms each being positive. This bound gets better 
with larger values of n, where kinetic contribution to Jia dominates, since each 'yai,ai+i is a 
better and better approximation to a geodesic on M'^. The injectivity radius is essential in 
reaching the approximation by M^. 

This gives the asymptotical growth of oscillator modes contribution to the winding string 
action, which in flat space we know is linear with respect to oscillator mode numbers. We 
can make this argument more explicit, for all n > we have asymptotic values 

/idK] = J2 /idK,,., J > J2 /id[7....,J >n-^- j = - (60) 

where the estimated energy are those from Dirichlet boundary problems with end points 
(o"j, o"j+i). Normalizing as usual the oscillator modes commutators, the energy of an excited 
state then grows linearly in number of oscillator operators. 

4.3 Diagrammatics (a further curvature expansion) 

The general calculation of the determinant requires knowledge of the metric on M. In the 
cases of interest, namely Ricci fiat manifolds, we are foxed to deal with the absence of such. 
We will here introduce a further expansion of the path integral exponent involving curvature 
tensors, which is a good approximation near the large volume limit. 

Order by order we still find terms which require the geometric information about M to 
evaluate, however, for a Ricci-fiat M we will see that the contribution linear in curvature 
vanishes. This lets us evaluate the path integral up to this order, with the consequence of 
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not including the local information ( 'germ') of the metric near the geodesic 7. This is a 
technical hurdle, that we hope to overcome in future work by studying approximately Ricci 
flat metricJ^. 

The leading few terms in this expansion are as follows 



fjp/lexp 



3 



dip' ■ dip' - Mjk{a) (^V^ 



Hp/ 



exp 



i=l 



(fa dip' ■ dip' 



1- / d:'aM,k{(T)ip'ip\a,t) 



(61) 



where we are working with free transverse bosons and treating the position dependent in- 
tegrated curvature term as perturbation. Inspecting the example of Eguchi-Hanson space, 
we find in that case the curvature term is explicitly diagonalizable over the whole geodesic 
(see I4.4p , if one could do so on a compact Ricci fiat space our result below can obviously be 
improved. It is not clear we can assume this scenario. 
In the free theory limit, we have correlation functions 



: ip'{a,t)ip^{a,t) :~ S'Hn\z - z\ 



(62) 



and so when M is Ricci fiat, the leading curvature expansion gives zero, since the trace of 
Mij vanishes. In such an approximation, we then find the NLSM torus partition function 
(here D = 4 but generalizations are obvious) 



^I-Ioop(t) 



E 

7 stable 



n 



(0;t; 



n 



=1,2 



=3 



(0;r) 



(63) 



for stable 7 all the frequencies are real and quantization of (free) oscillators is straightforward. 
The correction to this result will arise at curvature tensor squared ||-R^|| level, assuming the 
manifold M is Ricci fiat. This also discards interference between winding states, as can be 
seen from the simple sum over distinct 7's. When the closed geodesic has moduli, further 
modification is required. 

As can be seen from Appendix lA.ll we need to include additional interaction vertices 
in the worldsheet action to have a consistent curvature expansion. Examples are the cubic 
interaction for the transverse bosons and quartic interaction for the fermions. In a more 
complete calculation, at ||-R|P order we need to combine the term found at second order 



-"^^In case we had an explicit metric, we may simply quantize the string oscillators to be Q;Jj_|_^(cr)e*^""'"'^^*, 
<5Ji_i,(o')e''"~'^-'* where ol\+v^'^^ solves the curvature perturbed equation of motion. In effect, this is the same 
as reducing to one dimensional sigma model and evaluating the determinant of the index form. 
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here with the above mentioned terms which were dropped in our one-loop determinant 
approximation. These are comphcated polynomials in the curvature tensor, but they do 
have well established significance in the study of spectral geometry. 

There is an interesting connection between our expansion above and the so-called sub- 
principal wave-invariants introduced by [20] ^U\. The second order term in (l6Tll is an inte- 
gration of the worldsheet Green's function for free bosons. 

j (fz j £z [G'{z,z)fMjk{z)Mjkiz) = TT{M{z)G'{z,zfM{z)) (64) 

where the trace operation goes over both field indices and worldsheet coordinates. 

To make contact with the trace formula (for a particle), rewrite the Green's function as 
a sum over images 

G'{z, z) ~ — ^ \n {^\z — z + muji + nu2^) (65) 

m,n 

Now if we let the arguments in the Green's function go to zero and replace the divergent 
correlator by simply 5^ (a, a) and keep the sum over images. At | |-R| P order, after integrating 
over j d{a — 5"), we get an infinite sum with each term looking like the following 

J dt ^ da Mjk{(T)Mjk{(j) (66) 

This fits nicely into the residue expansion of the wave trace near the length of a periodic 
geodesic (i.e. its singular support). As explained in [10]. These are all local invariants 
depending only on the germ of the metric near the geodesic's tubular neighborhood and can 
be evaluated following the method surveyed in jlOj . 

The ^ da integrands of curvature polynomials performed in Fermi Normal Coordinates as 
we see above, are called Fermi- Jacobi-Floquet polynomials Ell in [ID]- Their integration 
along a closed geodesic orbit gives the 'subprincipal' wave invariants a^^k (see Theorem (5.1) 
in [inj)- The simple formula introduced in fl36p is the principal wave invariant a^^ associated 
to the trace formula. 

We won't give the explicit expressions for I^^k here, the explicit forms seem to be only 
known for two dimensional manifolds M (see eq. (5.5) in [lO]). The explicit construction 
of fc involves products of contraction of covariant derivatives of the Riemann curvature 
tensor with the tangent vectors ^ and with each other. This follows completely from a 
Fermi normal coordinate expansion, just as in our NLSM action expansion. 

In fact, since the functional leading to these wave-invariants are nothing but the length 
of 7 whose second variation is exactly what we call quadratic action for fiuctuations, we 
find it compelling to propose there is an exact agreement of the terms in our diagrammatic 
expansion and I^^k (when including higher degree vertex on the worldsheet). It will be 
an interesting problem work out the explicit forms of these I^^k curvature polynomials for 

^^The 'floquet' is simply the geometric series related to the Poincarc holonomy, explicitly defined as 
(3i = Y^-iW- for each angle of rotation. 
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general higher dimensional manifolds M. In terms of physics, the limit we are taking here 
is an optical limit of the path integral. 

4.4 An explicit example: Eguchi-Hanson space 

In D = 4 the obvious target space to consider is K?) which is Ricci flat. K?) has a moduli 
space which is locally a quarternionic coset, and the superconformal characters of K?) [26] 
provide an intriguing representation of the Mathieu group M24 [27tl28]. 

Constructing a geodesically complete exactly Ricci-fiat metric on K?) is obviously a chal- 
lenging problem on its own, which we will not discuss here. We will content ourselves at 
present with an explicit check of some of our basic ideas using an concrete non-compact 
Ricci fiat example. 

For this, we consider the Eguchi-Hanson space T*P^ [29] which provides a reasonable 
model of K?) [30] near each of its blow-up cycles, which may be situated in a family of 
Einstein metrics approaching the orbifold limit [31] . 

The geodesies on Eguchi-Hanson space are integrable, as shown in [32]. One would like 
to identify a stable geodesic on Eguchi-Hanson space. This is found by looking at the second 
variation, namely the index form (1551) for the chosen geodesic. As M is Ricci fiat, there are 
always negative sectional curvature components, and the stable geodesic if it does exist, is 
stabilized due to the kinetic contribution of the local basis around 7. 

Some simple cases are immediately clear, the blow-up cycle is geodesic with all its 
geodesies closed, namely all the large circles. A straightforward calculation shows that the 
curvature- mass term Mjj = Ryiyjl'-fla)) can be diagonalized leading to constant mass terms 
for the fluctuation modes. In particular, the longitudinal mode along 7 is indeed massless, 
there are two negative mass^ modes and one positive direction. The eigenvectors in the local 
frame rotate along 7, however the negative modes are not involved, so their kinetic energy is 
zero, and the large circles on are all unstable as expected. Additionally, the mass matrix 
entries have size and blows up in the orbifold limit which is a = 0. 

It is generally desirable to consider the possibility of a closed stable geodesic inside this 
small e-patch of the Kummer surface, we may look for orbits that either stays at a constant 
radius or oscillates between two turning points in the radial direction. The second choice 
can be ruled out in general for Eguchi-Hanson space. The argument is simple, and uses the 
explicit solution of the radial and ^-direction geodesic motion. We give here only the result 
of the calculations. 

Assuming that there is a oscillatory motion between two extremal values of r, then these 
must be points where f = 0. Indeed, the general case admits two such turning points. It 
also turns out that the 6 motion is completely determined by radial motion, as a result, 
we can find the change in angle 6 between these two turning points. We find that has 
non-vanishing imaginary part, which is unphysical. So all geodesies (coming in towards the 
bolt) must turn away at max(a, 2J/^/E). 

As a result, the only confined motion in Eguchi-Hanson space is a round trajectory at 
fixed radius tq = max(a, 2 J/v^)- Considering this case, we find the diagonalized mass 



24 



matrix to again have the same pattern of (0, — , — , +) eigenvalues. The nonzero values are 
{ — ^^1 "^4 ^) where has constant valuqlj. Again, calculating the combined value 

of kinetic (which vanish) and curvature terms, we find the negative directions stay negative. 
So these are also unstable geodesies. 

Another obvious choice is the radial motion. Due to the symmetry of Eguchi-Hanson 
space, we can do better, solving rather straightforwardly the case of motion in (r, 6')-plane of 
the polar coordinate. These are non-compact motions, which will escape to infinity. However, 
considered as part of the Kummer surface, they are interesting candidates for stable closed 
geodesies. Going to the inertial frame for the geodesic characterized by energy E and angular 
momenta J, the canonically normalized action for the fluctuating fields have three distinct 
radially varying mass^ values. 

They are the combinations of frame kinetic term and the original curvature induced mass. 
The former has a uniform value for all four eigenvectors of the curvature-mass matrix^ 

2, , AJ\r^-a^) 1 
^o(-) = - (67) 

This is the total mass^ for the eigenvector with zero Mij eigenvalue. This value is strictly 
positive, as r > a is always true. 

The two negative modes of Mij have identical values of mass^ 



4j2(r4 - a^) _ 2Ea^ 1 



= ""^^-^ -"^-ii (68) 



Clearly the first term dominates as large values of r, as the curvature has much faster decay 
as r~^. This showcases the general mechanism by which the geodesic stays stable against 
sectional curvature. If not for the first term, these directions are unstable directions and 
the geodesic would unwind along them. There is a region around the blow-up cycle which 
destabilizes geodesies coming too close to it. It would be interesting to study this further in 
a compact example, as what we see in this example could potentially work out there as well. 
The last mode has positive values for both terms 

"^2(r) = ^ + ^--, (69) 

For a crude model of compacitification we may cut a ball of radius A ^ a. The size of 
the integrated mass term in the worldsheet action is then 1/A^ <C 1 for the non-compact 
geodesic under discussion. As we scale down the Eguchi-Hanson space to fit inside a ball 
of size e cut out from say T^/1j2, the coordinates pick up inverse factors of e, making the 
most dominant term scale as e^^A^"'. This shows that so long we take the original ball to 



""^^When the orbit is at r = a, the energy is in fact E — and these values are infinite in the orbifold 
hmit as mentioned car her. 

^■^The radial coordinate r does not parameterize arc length of 7, but this is not significant for the current 
discussion. 
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be very large, the eventual integrated mass term, A ■ ~ e ^ 1 (or A^^ ^ 1), and can be 
safely neglected at leading order. The quantization of oscillator modes then resembles the 
free theory, leading to growth of density of oscillator states in the UV similar to fiat space. 

Summarizing, for the non-compact Ricci-fiat Eguchi-Hanson space, we find the un- 
bounded geodesies are stable against variations and the induced mass term agree with the 
picture that they do not cause significant change to the quantization of oscillator modes 
and their high energy level spacings. It is clearly desirable to test this picture further with 
a smooth closed geodesic on a compact Ricci-fiat metric on i^3, such as ones constructed 
using the gluing technique. 

4.5 Trace formula and modular invariance 

In this section, we will show that the trace formula (137|) transforms the worldsheet torus 
partition function ( l63l) in winding sector to that in the momentum sector. For simplicity of 
discussion, we take again dimension of M to be D = 4, while the general case is obvious. 

For the physics points we want to make, we want stress again that the trace formula is 
different from the modular property of the NLSM Hilbert space. The trace formula makes no 
assumptions about the stability of the geodesies included in the support of the wave-trace. 
In cases like the sphere, there are no stable geodesies to talk about, while the trace formula 
holds none-the-less. 

In the simpler case of Gutzwiller's trace formula [8], the 'empirical' distinction between 
stable and unstable CclSGS, clS also mentioned in [16], is that stable geodesies correspond to 
a 5-function peaks of the Laplacian spectrum distribution, while an unstable geodesic gives 
rises to smooth peaks with width. In this sense, while the trace formula does not discriminate 
against unstable geodesies, it also preserves its instability in the 'dual' spectrum. 

For a meaningful discussion of the Hilbert space of the NLSM, we need to talk about 
stable geodesies. In view of the 'preservation' of stability by the trace formula for a particle, 
we think of the trace formula as a machine that is useful for mapping true (i.e. stable) states 
in the winding sector to true states in the momentum sector. We will show that it indeed 
does this job, and comment on the evidence from this in support of our more general claims 
about the number of stable closed geodesies on a general Calabi-Yau manifold. 

First, we have the leading order partition function in the momentum sector, given by the 
Laplacian spectrum A = following from the trace formula. A genuine set of data for the 
(UV) asymptotics of the Laplacian spectrum on a Calabi-Yau is obviously more fitting for 
making a prediction for the geodesic spectrum, we do not have such data. In fact, numerical 
methods as developed in [53] [31] are more suitable for addressing the low lying eigenvalues 
of the Laplacian due to resolution limitations. This is the same limitation on being able to 
tell whether a discretized geodesic is closing (or closing smoothly). 

For a trace over the momentum sector of the Hilbert space, we will consider explicitly 
the bosonic contributions, as the fermionic case does not involve non-trivial zero modes and 
the corresponding discussion is simpler. The momentum states are given by a sum over all 
eigenvalues appearing in the formula (1371) and associated oscillators. 
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For the momentum sum, consider the (quasi-)eigenmodes corresponding to these asymp- 
totic eigenvalues are Gaussian beams, which is of the form 



^a(X) ~ e 



iXa+eicT)MijicT)ma) 



co(0 H- H- 



(70) 



and can be inductively solved (see [IH]), with A the eigenvalue and Mij{a) the curvature in 
Fermi normal coordinate. Vertex operators for momentum states can be constructed from 
the eigenf unctions of Laplacian by adding oscillators (derivatives of X), for example 

V'^iX)^dede^xiX{a,t)). (71) 

To leading order (where the left and right-moving Hilbert spaces decouple), we can take 
the bosonic partition function in the momentum sector to be 



Z_(^)^tr(.V°)^ tr(.V)l.(r)r ^ 



nil 
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(0;r) 



(72) 



where 



A = (2?™ + aimi + a2m2 + a^m^ + z/^)^ + o{n ^/^) 

(L, 



^1) 



A = ^ {2-1171 - aiTTii - a2m2 - asm^ - u^f + o{n ^/^) (73) 

where we have renamed rotation angles to avoid confusions. As before, we have ra, mj G Z, 
and q = e^'^*^. Notice that by choosing the form of A we are specifying the spin of the 
momentum states. 

Notice that invariance of Zmome.i^) under r — )■ r + 1 is not obvious for generic values 
of the rotation angles which are irrational. This can be remedied by either making the 
momentum sector vacuum picking up a phase proportional to Ti or introducing an explicit 
Wilson line for the rotating frame. In the former case, the vacuum is aligned to the creation 
and annihilation operators which rotates around 7 (see 13. 3p . and Hermiticity of the action 
is preserved. 

Next we show that this sector of the partition function over momentum states transforms 
under r — )■ — ^ into the partition function for the winding states. We apply the Poisson 
resummation formula to the sum over n, this turns the momentum factor into 



tr(gV) = ^J2Y1 (e-^''+^("-+^' e'^^'"'^^'"''"'^'^ (74) 



Recall the familiar transformation T2 ^ and let us focus on the first term in the exponent, 
then 
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giving precisely the classical contributions to the partition function of iterates of a winding 
string of length L^. 

The other factors in ( !74l) deserve a separate explanation, they transform under r — )■ — l/r 
into 

e^2^ T^^j e ^ (76) 

The second term arises from the shift of zero point energy by twisted boundary conditions 
as in (145|) . with holonomy induced shift by {y^ + ^-mjai). When we include fermions, 
these terms cancel in the Neveu-Schwarz sector. The ^ term comes from the spin of the 
momentum sector states, where the r — )■ r + 1 invariance was restored by frame rotation 
under P^. 

To compare the rest of the momentum sector partition function with the winding partition 
function, recall for transverse bosons with nontrivial holonomy in the winding direction 
(spatial twist) 
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(0;r) 
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And the total bosonic contribution to the winding sector partition function from (163] 



(77) 



Z wind.{T) — 



nil 
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(7J 



Combining (172|) . ( fr4|) . ( !76|) and (178|) . and up to factors of y |t| which can be easily fixed 
from the number of degrees of freedom, we see that the nontrivial part of the transformation 
under r — )■ — 1/t exchanges the temporal and spatial twists exactly as follows from Z^^ij) = 
Zl^_a{—l/r). It's also useful to realize that complex conjugation reverses the sign of the 
characterestics in the ^-function. The fermionic case follows exactly the same pattern, with 
proper modifications of periodicity for the ^-functions.. 

We can conclude from the above calculation that the trace formula, when applied to 
stable closed geodesic spectrum, or equivalently the associated Laplacian eigenvalues in the 
sense of [TB], does lead to modular invariance of the torus partition function in the limits 
under consideration. This however will not prove our general claims as put forward in the 
introduction, which has additional motivations from well known facts about conformal field 
theory associated to Calabi-Yau spaces and their orbifold limits. The r — >■ — ^ transform of 
the momentum sector partition function which itself is more straightforward to obtain, is a 
more stringent check of the limit we used to obtain the winding sector partition function. 

Still, it is a mild assumption that the Laplacian spectrum gives the momentum sector 
of the NLSM (at least in large volume limit), and a further reasonable assumption that the 
closed geodesic contribute as winding states to the NLSM partition function. Then based on 
the asymptotics of the Laplacian spectrum, which can be obtained entirely independent of 
the trace formulas, and has definitely polynomial growth, it seems the stable closed geodesies 
on a Calabi-Yau manifold should have the same corresponding polynomial growth as L^. 
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From the mathematics hterature, the discussion of stable closed geodesic is still a chal- 
lenging problem. There is the famous result of Serre [37] that on a compact Riemannian 
manifold there are infinitely many geodesies connecting any pair of points. The proof is 
based on spectral sequence techniques and makes essential use of the path/loop space AM 
(and QM). Later, Gromov's result [38] give lower bounds relying on information of homology 
of the path space, the number of geodesies between a pair of points p,q E M is bounded 



where a and b are (which are not algorithmically given) constants and 6„(AM) the betti 
numbers. For dominantly most manifolds, the above sum grows exponentially with respect 
to L (examples are discussed in jllj). 

In rare cases, the growth with respect to L is polynomial [39] and this includes elliptic 
surfaces such as K3 (we are not aware of claims about Calabi-Yau 3- folds). In fact the growth 
applies also to closed geodesies, in the cases mentioned. This is circumstantial evidence for 
our proposals, however there are not sharp statements (aside from ours) about the power of 
the polynomial growth with respect to length. 

In view of this, further mathematical evidence, presumably from studying approximately 
Ricci flat metrics asymptotic Laplacian spectrum, would form very interesting inputs for 
establishing (or disproving) the asymptotic behavior of geodesic length spectrum we propose. 
Finding the explicit (2nd variation-)stable closed geodesies will of course add substantial 
support to our present proposals. 

Incidentally, there is a well-known Splitting-theorem [40] prohibiting the existence of 
globally length-minimizing geodesies on manifolds which has nonnegative Ricci-curvature 
(including Ricci-flat) and is not of the direct product form M^-i xM. The locally length- 
minimizing geodesies we study, which are stable due to positive second variations of the 
action, are not banned by this theorem. 

5 Comparison with orbifold CFT 

In this section we make some exploratory comparison with the orbifold limits of K3. Our 
main goal is limited to showing the consistency of the (free) CFT partition function with the 
presence or absence of winding sectors in the large volume limit, on Eguchi-Hanson space. 

We consider the non-compact orbifold M'^/Z2 which is the singular limit of the Eguchi- 
Hanson space discussed in (14. 4p . The discussion is nearly identical for T^/Z2. 

The action of Z2 simply inverts all (signs) of coordinates, the action on the (free) 
fields are then 




(79) 



n<bL 



(80) 



At least to mathematicians and we thank Michael Anderson for clarifying this point for us. 
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The action on fermions follow from preserving M > (2, 2) worldsheet left and right moving 
supersymmetries. There is an unique fixed point (the origin) for ]R^/Z2 while there are 16 
for TVZa. 

For relation with the winding states, we first consider the untwisted sector on the orb- 
ifolds. Here we have indeed states which carry momentum, identical to the double cover 
or torus. The untwisted sector has the partition function 



■^untwisted (''") 



1 



^(0,0) 



^(0,1) 



V)] 



(81) 



where 



^(0,0) 



(r) 




p2 -y 



4 2 



9-Ar) 



9,{r) 



t]{t) 



7]{t) 



(82) 



is modular invariant. This is the only sector containing zero modes (momentum) which are 
projected out in the twisted sectors. 

In case of T^/Z2 the integration over the continuum is replaced with a sum over momen- 
tum and winding, as in the tori case, and as in that case modular invariance is manifested 
as the r — )■ — 1/r exchange of momentum and winding sectors. These we expect to turn into 
stable closed geodesic once we turn or marginal deformations, which is in the twisted sector. 
Notice that the set of closed geodesies on tori has measure zero. 

Modular invariance of the second term of ^untwisted (with insertion of Z2 element) dictates 
inclusion of twisted sectors, and the following combination is required 



^(0,1) (^) + ^twisted = ^(0,1) (t") + ^(1,0) (^) + ^(1,1) (^) 



(83) 



where the subscript (r, s) is the usual notation for twisting in time and space direction of 
the worldsheet. Each sector has a 16-fold degeneracy which we divide by and a weight of 8 
from the orbit, so we find an overall factor of | as in ( IHTj) . 

For the (r, s)-twist sector, consider the bosonic and NS-sector fermion contributions (for 
J\f = (4,4) one can include the SU{2)i isospin) 



^3 {Z + ^) ^3 {Z - ^) 



61 (^) 

The R-sector result comes from inserting (— )^ and shifting z ^ z + ^i^. 



from [24] one can rewrite the holomorphic factors 

'9^iz)Y (e,{z] 



Xr 



0; 



9, 



' r+sr ^ 



(84) 

Using formula 
(85) 



Each of these contributions can be decomposed into the sum of a BPS character and 
infinitely many non-BPS (massive) characters. Explicit forms of these decompositions can 
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be found in e.g. [23]. For the untwisted Z(^o,i) the BPS character is that of the identity 
representation with {h, I) = (h, I) = (0, 0) . For the twisted sectors, the BPS representation 
(in NS sector) is that of {h,l) = {h,l) = (1/2, 1/2) which is marginal and corresponds to a 
moduh. In the c = 6 theory with Af = (4, 4) worldsheet superconformal algebra, these are 
the only two possible BPS representations. 

The twisted states are realized by the vertex operator involving the unique twist field 
a{z, z) of dimension (1/4, 1/4) in the case of M^/Z2 

and for T^/Z2 one has 16 such twist fields. The twisted partition function in NS sector has 
the decomposition [35] 

^twisted (r, y) = Ixi.oP + 1X1,1^ (87) 
y = e^'^*^ and explicitly the holomorphic blocks are 



Xi,o 



^3 ; V xe^ir 

-1/2 



03 J \ 03 J ye-sir) 

2 //i/N\2//i/x\2 



On the other hand the BPS and massive characters are given by 



73 

r/(r)3 

mi 
q 2 8 



r/(r) 



hir) 



V{r)03{rj—i + q 



~1 1 + 2 



where h^lr) is one of the so called Mordell functions. To expand the holomorphic blocks 
in terms of the characters, we use the well known identity 2rf{T) = 02{t)63{t)64{t) , which 
leads to 

X,,o = cK^ (1 = 1/2) -(hir) ' W^^^3(.)^^ 



Mr)y \vir) 

Xi,i = ch^^(/ = l/2)-/i3(r) 



7]{t) 

. ehr(/ ^ 1/2) - (/.3(r) H- ,4Mlll) (?My (80) 
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Expanding in q the coefficient of then leads to an infinite sum over massive char- 
acters. All three series can be considered a power series of and share the lowest order 
term of gs, which gives the lowest operator's dimension = i. For more details on these 
expansions, see for example 

It is curious the massive characters q -^^^ appeared in our one-loop determinant in (1521) 
when the massive transverse bosons are decoupled. This is probably just a coincidence. Can 
these twisted sector massive representations can be associated to winding strings, on K3 or 
Eguchi-Hanson spaces? Due to Z2 projection, in a twisted sector the winding string must 
lie in a fixed locus, for example the interval connecting two fixed points. 

Based on the fact that these characters appear for the non-compact ALE-spaces, espe- 
cially Eguchi-Hanson [35], this possibility seems unlikely, since for the latter space there is 
only one unique fixed point. Furthermore, in the semi classical calculation leading to (1521) , 
to really decouple the transverse bosons the curvature must be large along the whole length 
of the winding string, which isn't realized unless the winding states are on the blow up cycle 
that shrinks in the orbifold limit, which we know as unstable. 

As far as Eguchi-Hanson space is concerned, we find it consistent with the known pic- 
ture where twisted sector states are associated with the Kahler modulus and the massive 
contributions as observed in [35] simply complete the BPS characters r(2) invariant. 

One (of many) question that seems to arise from these considerations, that of explaining 
how the rather transcendental spectrum of the closed geodesies could evolve into the rather 
rational operator spectrum of an orbifold CFT, is exemplified in the above consideration 
of massive characters for the twisted sector. This was one of our original motivations for 
undertaking the studies in this paper and it remains a mysterious question, while we hope 
we've made clear the direction we believe could lead to its answer. 
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A Appendix 

A.l Fermi normal coordinate expansion 

We start with the bosonic action, which reads 




(91) 
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usually T = Xjli^ol . hai3 is the worldsheet metric and G^^ that for the target space. Here 
indices are = cri,(J2', and fi,^ are space-time coordinate labels. We find the harmonic 
map equation in explicit form 

Vh-^da{Vhh"^dpX^) + h''^Tl^^d^X''dpX^ = (92) 

here F^;^ is the Christoffel symbol for the metric G/^jy. In conformal gauge, the equation 
reduces and for solutions depending on a only we find the geodesic equation. 

We make use of Fermi normal coordinate and expand around a closed geodesic 7(0") in 
target spacj^. For complementary details see |41]. The Fermi conditions are 



9fiuU{t) = Sf,u , ^uxUit) = (93) 

Taking a geodesic Ao-(s) going off the closed loop 7(0") parametrized by s, we set up a 
coordinate chart on a open set. Using the solutions of the geodesic equation we can Taylor 
expand the worldsheet fields 

As y^{s) is a geodesic, applying the geodesic equation leads to 

x'^is) = x^{o) + c^s - ^Kx\s=oCes' - ^ (rf,,,,) - 2rf^,rl^)) \s=oCees' + 0{s') (95) 

Using ( p3|l one can show also the completely symmetrized partial derivatives of the Christoffel 
symbol in the transverse directions vanish, for example 

ri;M)l7 = (96) 

In the Fermi normal coordinate 

X\s) = X\0) + Cs , X°(s) = a (97) 
Taylor expansion once again gives the equation for the local basis 

Taking the initial basis to be orthonormal e^^^(O) = 6'^^ and solve the transport equation 

de'^ 

Del)ldt=^ + Tl^el^u^ = Q (98) 
We recover the Riemann normal expansion 

(^) = + ^^%oo f + ^i?%oo;0 (99) 
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FNC also exists for accelerated or rotating observers. 



33 



Take the base point for geodesic line A(s) on 7(0-) to be identified as s = 0, the geodesic 
equation solution gives the coordinates of geodesic motion 



d 1 d 

2 
1 



= C''W+2^"ei;|Ws+gfl%o(« = 0)(x'ef„(iyej',,(i))s'i+... (100) 

Notice here the initial velocity is given as x*e|'^^(t) in the basis determined by f l99|) . addition- 
ally the Christoffel symbol r^;j^(t) is expanded similarly along 7(t). Expanding out to first 
two orders, we see the following results 

C%s)=t + ^RoiOj{x'x^)sH + ... 

C{s) = t + iK.oi^nst' + ii?},o(a:V>2t. . . (101) 

The fermi normal coordinates for the above point is = (^j^;*), so the coordinate 

change required to fermi normal coordinate is 

c^s) = e + lRo^oAee){e) + ... 

c{s) = t + lRi>,o{mr + iRU^^x'^m^em . . . im 

And applying the transformation rules (y'^y(FNC) = gxp{KNC)^^ we finally arrive at the 
expansion of the metric in fermi normal coordinate 

goo = i+ Ro^oJee + o{e) 
goi = R^koi^e + o{e) 

g,, = 5,, + Rik.ii^e + 0{e) (103) 

The bosonic nonlinear sigma model action is expanded order by order in Fermi normal 
coordinate, where is the background field much like the expansion in Riemann normal 
coordinate in |12]. Notice at fourth order in we have terms involving Riemann tensor 
squared. Using (j97|) the action reads 

/[X + tt] = j d'z^S^^d^X^d^X'' 

+ J d\ i{(5,,v,ev,e + (v«eo)' + i?o^jo(t)9aXUx°ee} 

+ j d^z{^Ro,Ut)^ae + Rojkom»^o)dc.X''^^e 

+ J d^z ,(t)9,x09„xoeee' + • • • (104) 
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Conformal condition easily follows. Wave function renormalization is identical for the 
background and fluctuation fields. The only non-zero interaction term to this order is a 
spacetime dependent mass term 

Roijoit)d^X''d^X'ee 
Contribution to the stress tensor is 

a'RoiojS^^ ln(A//i) = a'Roo ln(A//i) 
The fixed point requires the component Rqq of Ricci tensor vanishing. 

A. 2 Fermion in Fermi normal coordinate 

The two dimensional JV = 1 supersymmetric nonlinear sigma model has the following action 
in components form 

I[X,ij] = i y d'zg^,{X)d^X>^d^X'' + tg,,{X)^^rD^r + ^R,uXpiVr)i^'r) (105) 

We use the following notation for the Fermionic fields and gamma matrices in two Eu- 
clidean dimensions 

f ^ 1 

where C is the charge conjugation matrix. The action then reads 

I[X,ij] = d^zg^,{X)d,X^d,X'' + lg^,{X)iJlD,r+ + ^9AXW-D.r- + \R^^^^^ 

(107) 

where D^^^^ = 9^^^^ + r(;^(X)a„X>^. 

In the Fermi normal coordinate, we have the following Taylor expansion of the Christoffel 
symbols along the geodesic 7(t), where x° = cr, = ^*s. 

The worldsheet supersymmetry transformations are 

SX^" = eilj" 

6^;^ = -$X^e-r''^^{e^'')^^ (108) 



= (? J)' ^' = "^=(o -l) 

(° (106) 
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In components they read 

= -e^d,X^-ie_r^Yl^^^^ (109) 

When the target space manifold is Kahler, there is an additional supersymmetry invari- 
ance under 

6wr = -$x^^^ + h:^^^j%j\{rry (no) 

The complex structure J^p allows to split the complexified tangent bundle T^X = T^'^X © 
T^'^X, and the extended A/" = (2, 2) supersymmetry transformations parameterized by the 
pair of spinors (e, e) read in complex coordinates 



6X' = 




SX' = 




6iji = 




Si^l = 








= 





(111) 

We shall note here that in supersymmetric theories for localization arguments to work, one 
needs the Fermionic fields to satisfy periodic boundary conditions since the supersymmetry 
transformation parameters do. 

Collect the classical and quantum terms of the action. Bosonic terms 

4 = / d\U,,,d^X^d^X- 

+ y i{5,,v,ev„e + (v„eo)' + i?o^,o(t)5a^°9a^°ee} 

(112) 

Fermionic terms are important only to quadratic level 

Il = ^jd'zz {^iD.^l + ^° ) + t {iP\D,^l,\ + ^ID^^l) (113) 

Higher order terms in the connection give rise to interaction vertices between bosons and 
fermions and are dropped. 



36 



A. 3 Geodesies break worldsheet SUSY 



First we use the real notation for M = (1,1), where it makes sense to separate the parallel and 
transverse fields. Along the geodesic, the affine connection vanishes in the normal coordinate 
and the supersymmetry transformation rules are 

= -e+d.X^" (114) 

In the fermi normal coordinate, the geodesic is simple and the worldsheet can wind integer 
times around the geodesic length L^. Let the periods of the worldsheet torus be ui and U2, 
the generate the lattice A such that S = C/A. Then denote the modulus by r = U2/u}i and 
the boundary conditions for the zero-modes are (without Poincare map effect^) 

X\ai + Re{ui), 02 + lm{bJx)) = X'{ai + ReM, (J2 + Im{uj2)) = X'{(Ti, aa) 
X°(ai + Re{uJi), ^2 + Im{uji)) = X°((Ti, ag) + riiL 
X°(ai + Re{u2), (J2 + Im{u2)) = X°((Ti, aa) + ngL 
tp^{ai + Re{uji), 02 + lm{bJx)) = ±tp''{ai, aa) 

^^(ai + ReM, a2 + Im{uj2)) = ±V{(^i, (^2) (115) 

The zero mode part of the worldsheet fields are determined from the above boundary 
conditions, especially 

XOf(j J. \ _ nilm{uj2) - n2lm{uji) niRe{u2) - n2Re{ui) ^ 

Re{ui)Im{u2) — Re{u2)Im{ui) Re{ui)Im{u2) — Re{u2)Im{ui) 

(116) 

Recall our convention that z = (cti + ia2)/2, z = (cti — ia2)/2 this gives 

X%z, -z) = -"^f L . - . "f^-"^"^ L .- + ... (117) 

Im[uJiUj2) Im{UiU2) 

The ground state associated to the classical geodesic solution has the following non-zero 
variations when acted on by supersymmetry 

Si/Jl = -e-d^X" = -I ^Le- 

Im[UJiUJ2) 

. ,0 o .niU2-n2UJi 

dip'i = -e+d^X" = I ^Le+ (118) 

im[ujiuj2) 



^^The discussion here only concerns the classical solution. Further, as we argue in the main text, the 
holonomy around the closed geodesic does not effect the super-current. 
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Here e± are the two independent real supersymmetry transformation parameters. Trading 
e± for the pair of hnear combinations (ei, €2) = C^^"^*^' , ^^~''' ), we have 

Im[UJiUJ2) 

. ,0 , •'^1^2 - n2U}i ^ 

Im[uiiui2) 

,0 .niuj2 - n2UJi 
diip- = I J ^i^ei 

Lm\L>JxL>J2) 

02W- = I J , ^^£2 (119) 

Lm\L>J\L>J2) 

Obviously one linear combination of the pair ('?/''^, ) is invariant under transformation 
parameterized by e\ and another linearly independent combination is invariant under €2. 
Explicitly we have 

bx{{nx'^2 - 'n2'^^'\\)\ + {nxL>j2 - n2io^iPj\ = 

62[{niUJ2 - n2UJi)ip1 - {niU2 - n20i)ilj'^] = (120) 



We see that generally they are independent, except when one of the coefficients vanish which 
corresponds to dzX^ = or dzX^ = 0, i.e. anti-holomorphic and holomorphic maps, which 
leaves respectively ip^ and invariant. 
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